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ABSTRACT 

We place constraints on the average density (f2 m ) and clustering amplitude (erg) of matter using a combi- 
nation of two measurements from the Sloan Digital Sky Survey: the galaxy two-point correlation function, 
w p (r p ), and the mass-to-galaxy-number ratio within galaxy clusters, M/N, analogous to cluster M/L ratios. 
Our w p (r p ) measurements are obtained from DR7 while the sample of clusters is the maxBCG sample, with 
cluster masses derived from weak gravitational lensing. We construct non-linear galaxy bias models using the 
Halo Occupation Distribution (HOD) to fit both w p (r p ) and M/N for different cosmological parameters. HOD 
models that match the same two-point clustering predict different numbers of galaxies in massive halos when 
f2,„ or a% is varied, thereby breaking the degeneracy between cosmology and bias. We demonstrate that this 
technique yields constraints that are consistent and competitive with current results from cluster abundance 
studies, even though this technique does not use abundance information. Using w p {r p ) and M/N alone, we find 
fi£ 5 cr 8 = 0.465 ± 0.026, with individual constraints of O m = 0.29 ± 0.03 and <7 8 = 0.85 ± 0.06. Combined with 
current CMB data, these constraints are fi m = 0.290 ±0.016 and cr 8 = 0.826 ±0.020. All errors are la. The 
systematic uncertainties that the M/N technique are most sensitive to are the amplitude of the bias function of 
dark matter halos and the possibility of redshift evolution between the SDSS Main sample and the maxBCG 
cluster sample. Our derived constraints are insensitive to the current level of uncertainties in the halo mass func- 
tion and in the mass-richness relation of clusters and its scatter, making the M/N technique complementary to 
cluster abundances as a method for constraining cosmology with future galaxy surveys. 

Subject headings: cosmology: observations — galaxiesxlustering — galaxies xlusters 



1. INTRODUCTION 

Galaxy clustering measurements offer a unique window 
into the distribution of dark matter in the univers e. The re- 
cently completed Sloan Digit al Sky Survey (SDSS: lYbrk et al.l 
l2000l lAbazaiian et al.l 120091) provides unprecedented preci- 
sion and accuracy in its map of the local universe. However, 
one of the main impediments to the use of galaxy cluster- 
ing for inferring cosmological information is galaxy bias — 
it is probable that the distribution of galaxies differs from 
the distribution of dark matter. Thus both the amplitude and 
shape of the galaxy clustering signal are biased relative to 
the clustering of dark matter at quasi-linear and nonlinear 
scales. This bias is degenerate with cosmology in such a 
way that a bias model can be constructed to match the ob- 
served real-space galaxy two-p oint correlation function for a 
range of cosmological models (|Yang et ai1|2004 ; iTinker et al.l 
l200l IY00 et al.ll2006l: IZheng & Weinbergfl2007l) . In this pa- 
per, we combine real-space galaxy clustering and the cross- 
correlation of galaxies with galaxy clusters. The mass scale of 
the clusters is determi ned by weak gravitational lensing from 
ISheldon et al.1 (l2009bl) . which provides the leverage to break 
the degeneracy and constrain cosmological parameters. 

The framework within which we construct the bias 
model is the Halo Occupation Distribution (HOD; see, 



e.g., Peacock & Smith 2000; Seliak 2000; Scoccimarro et al. 
200 It iBerlind & Weinberg! 120021: ICooray & ShethlL2002l In 
this model, all galaxies are contained within dark matter 
halos with a probability distribution P(N\M), the probabil- 
ity that a halo of mass M constrains N galaxies of a given 
class. We define a dark matter halo as a collapsed, virial- 
ized object with a mean interior density of 200 times the crit- 
ical density. The statistics governing the spatial distribution 
of dark matter can be estimated with analytic models (e.g.. 
Press & Schechterl [19741; ICole & Kaiserl 119891 iMo & White 



1996 t ISheth & Tormenl 11999b ISheth et al.l l2001t lMa"et"a! 
i2010t> . but they are most accurately calibrated through N- 



:curately 

body simulat i ons dJenkins et ailllOOUfSeliak & Warrenll2004 
Warren et al.l l2006t iReed et al.l l2007t ITinker et al.l 12008a; 



Pillepich et al Miarnriker' et al.ll2010T) . These simulations 



specify the mass function of dark matter halos, the large-scale 
bias of halos, and their quasi-linear clustering relative to the 
clustering of the matter distribution itself. Thus, once the 
HOD is constructed with a model for the internal distribution 
of galaxies within halos, this relationship between galaxies 
and halos fully specifies the spatial distribution of galaxies 
on all scales. The key ingredient within P(N\M) is the mean 
number of galaxies as a function of halo mass, (N)m- Dis- 
tinct c osmologies produc e distinct populations of dark matter 
halos (IZheng et al.ll2002l) . thus for each cosmology there is a 
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distinct P(N\M) and (N)m such that the model matches the 
same galaxy two-point correlation function £ g (r). 

Data sets that are sensitive to the underlying mass scale 
of dar k matter halos can break the bias-cosmology degen- 
eracy (IZheng & Weinberg! 120071) . The most direct measure- 
ment is the mean number of galaxies in halos of mass M. 
In this paper we will use the ratio of these two quantities, 
M/N. This measurement is most easily made in the most 
massive halos — galaxy clusters. Cluster-sized halos are the 
largest collapsed structures in the dark matter density distri- 
bution, making them relatively easy to locate observationally 
and providing myriad methods with which to estimate their 
masses. The M/N measurement is analogous to the mass-to- 
light ratio of dark matter halos, M/L. In previous decades, 
M/L of galaxy clusters was utilized as a met hod for inferring 
the dark matter density parameter £l m (e.g.. | Gott et al.l | 1 974 ; 
Peebleslll986tlBahcall et al1fl99l 12000: ICarlberg et al.|[l996 : 
Rines et al i EoolT The pivotal assumption in this method was 
that the mean M/L of clusters was representative of the mean 
M/L of the universe, thus fl m = (M/L) x p\ um / /? cr it. where p\ am 
is the luminosity density and p CI i t is the critical density. These 
results were among the first to challenge the prevailing theo- 
retical expectation that Q, m = 1 . However, efforts to use this 
relation yielded values of il m near 0.1, in significant tension 
with a large array of other methods that were converging on 
a 'concordance' cosmology of fl,„ = 0.3 and as = 0.9 (e.g ., 
ISpergelet al.ll200l iTegmark et all 120041: ISeljak et al.ll2005l) . 
However appealing, the assumption that clusters are repr esen- 
tative had little theoretical support (Ba rdeen et al.lll986h . and 
the lack of agreement between cluster M/L and other mea- 
sure s of f2 m strongly imp lied that they are biased objects. 

In iTinker et al.l d2005l) . we demonstrated that cluster M/L 
ratios depend not only on the matter density but also on the 
amplitude of density fluctuations, erg, once the galaxy bias 
model is constrained to match the projected galaxy corre- 
lation function w p {r p ). Combin ing M/L measureme nts of 
clusters in the CNOC2 survey (ICarlberg et al.l 119961) with 
early clustering results from SDSS (IZehavi et al.l 120051 ). we 
determined that these data constrain a degeneracy curve 
(fi m /0.3)°' 6 (cr 8 /0.9) = 0.75 ±0.06. This degeneracy is sim- 
ilar to that constrained by measurements of the abundances 
of galaxy clusters; 'cluster normalized' models that lie 
on the same value of fibers predict the same number of 
cluste rs, where 7 is usually around 0.4 - 0.6 dRozo et ail 
120101). These results were i n good agreement with those 
of Ivan den Bosch et alj ( 120031) . who combine d clustering in 
the Tw o Degree Field Galaxy Redshift Survey dNorberg et al.l 
120011) w ith cluster M/L and first-year CMB resu lts from 
WMAP dSpergel et al.l2003l) . IVale & Ostrikerl (120061) reached 
similar conclusions using a subhalo abundance matching 
method, again drawing on cluster M/L ratios. All of 
these studies demonstrated that the 'concordance' values of 
(fibers) = (0.3,0.9) could not simultaneously account for the 
clustering of galaxies and the M/L of clusters. These conclu- 
sions were validated when three-year WMAP re sults emerged 
in ex cellent agreement in the Vl m -a% plane dSpergel et al.l 
120071) with the M/L resu lt offfl ,„/0.3) x (cr 8 /0.9) b - 6 = 0.75 ± 
0.06 from ITinker et al.l (120051) . The key revision to the 
WMAP estimates was improved correction of polarization 
foregrounds, which reduced the estimated electron scattering 
optical depth, and which thereby reduced the matter fluctua- 
tion amplitude inferred from the CMB anisotropy amplitude. 

Although the confirmation with later CMB data 
demonstrated the potential of this approach, both the 



Ivan den Bosch et all d2003l) and ITinker et al.1 (120051) studies 
paid little attention to the possible systematic errors in 
the use of halo occupation methods, emphasizing instead 
the systematic uncertainties in cluster M/L estimates, this 
paper, we use the largest data sets currently available. For 
the clustering we use measurements from Data Release 7 
of the SDSS dAbazaiian et alj 120091) . which covers nearly 
a quarter of the sky to obtain close to a million galaxy 
redshifts. Our cluster sample com es from the maxBCG 
cluster catalog (Koester et al. 2007a), consisting of over 10 4 
optically-detected clusters. These massive data sets close 
the door on statistics-limited analyses; in this paper we take 
great care to quantify the systematic uncertainties in our 
theoretical model. Our knowledge of the statistics of dark 
matter halos is not perfect, and our assumptions for how 
galaxies populate halos are not infallible. Thus we will 
incorporate these uncertainties into our analysis in order to 
yield robust constraints on cosmological parameters. 

This paper is structured as follows: In section 2 we present 
our clustering data and our measurements of M/N from the 
maxBCG clusters. In section 3 we describe our theoretical 
model for both w p (r p ) and M/N within the HOD framework, 
demonstrating how this method is sensitive to different cos- 
mologies. In section 4 we present a detailed description of 
the systematic uncertainties in our analysis. In section 5 we 
present our results, showing constraints in the £l m -cr% plane 
and constraints on HOD parameters. In section 6 we discuss 
these results in context with other studies that utilize clusters 
as a cosmological tool and present prospects for the future 
with the M/N method. We briefly summarize in section 7. 

2. DATA 

2.1. Galaxy Clustering Measurements 

Our galaxy clustering measurements are obtained from 
7300 deg 2 of sky, slightly smaller th an the full SDSS DR7 
of 7900 deg 2 dAbazaiian et alj|2009l) . The methodology is 
described in detail by IZehavi et al.l ("2010: hereafter Z10). 
Briefly, we choose three volume-limited samples of galaxies 
on which to focus our analysis: 01 M r < -19.5, 01 M r < -20.5, 
and 01 M r < -21.0 (also listed in Table l) 1 . These sam- 
ples have complementary attributes for this analysis. The 
M r < -20.5 sample, which corresponds approximately to a 
sample of all galaxies brighter than L», is the optimal sample 
for this type of analysis; the magnitude limit is bright enough 
such that the volume contained suppresses sample variance in 
the clustering measurement, but it has a high enough number 
density that the (N) m will be above unity for the mass range 
probed in the sample of clusters we use. The fainter sam- 
ple, M r < -19.5 maximizes the signal-to-noise in the M/N 
measurement, while the brighter sample, M r < -21.0, has 
the smallest sample variance for the clustering measurement. 
The galaxy clustering measurements are only weakly corre- 
lated because z max is different for each sample, minimizing 
the overlapping volume (see Table 1). 

The clustering quantity we utilize for each sample is the 
projected two-point galaxy autocorrelation function, w p (r p ). 
This statistic is defined as 

1 We do not assume a value of h for galaxy magnitudes, but for brevity we 
write absolute magnitudes as M rather than M — 5 log h. The superscript indi- 
ca tes the redshift to which t he magnitude is ^-corrected using the algorithm 
of lBlanton & Roweisl UMfo . 
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FIG. 1 . — Panel (a): The background-subtracted projected galaxy number density around maxBCG clusters of varying richness (see key in panel [b]). These 
measurements are for galaxies brighter than 25 M, = —19.0. Panel (b): The deprojected number density profiles of galaxies for the data in panel (a). The filled 
circle indicates the halo radius Rn200 for each richness bin. Panel (c): The variation in « gii i(r) as a function of magnitude threshold. Panel (d): The M/N 
measurements as a function of cluster richness. The data are labeled by their °'M r magnitude thresholds (see Table 1 for the conversion between 0A M r and 
25 M r ). Note that we define N to include on satellite galaxies. These data do not include the correction for miscentering, discussed in §4. 



w p (r p ) = 2 / ^(r pi n)dn, (1) 
Jo 

where r p is the projected separation between two galaxies, 
7r is the line-of-sight separation between two galaxies, and 
£(r p ,ir) is the measured two-dimensional correlation function. 
Due to the peculiar motions of galaxies against the Hubble 
flow, £(r p ,7r) is anisotropic. Integrating along the ir direction 
minimizes the effects of redshift space anisotropics, allowing 
for an easier comparison to analytic models. If 7r max = oo, 
redshift space effects are eliminated entirely, but the finite 
volume of the survey requires that we limit the integral to 
TTmax = 40 h~ l Mpc. At this 7r max , redshift-space effects on 
w p {r p ) are larger than the error bars on the DR7 measure- 
ments at scales of r p > 5 /r'Mpc. Although it is possible 
to include redshift-space anisotropics in HOD calculations 



of the non-linear two-point correlation function (e.g., iTinkeri 
l2007h . we limit our analysis to r p < 3 h 1 Mpc to eliminate 
any systematic biases due to redshift- space effects . Our the- 
oretical model for w p (r p ) utilizes the ISmith et all (120031) fit- 
ting function for the nonlinear matter £(>), thus limiting our 
analysis t o small scales elim inates biases from possible er- 
rors in the ISmith et ail (12003b prescription. Measurements at 
these scales provide more than enough statistical power for 
our constraints given the systematic uncertainties discussed 
in §4, which dominate our final error budget. 

We use the full covariance matrix for each clustering sam- 
ple. The covariance matrix is estimated through the data us- 
ing the jackknife technique. The full survey is divided into 
104 subsamples of roughly equal sky area (thus nearly equal 
volume), and each subsample is removed from the survey and 
w p (r p ) is re-measured. The jackknife cumulative covariance 
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FIG. 2. — The evolution of the r-band cumulative luminosity function 
between z = . 1 and z = 0.25. The black curve shows the measurements from 
IBlanton et al.l 120031) . The dashed curve rep resents our estim ate of the °' 25 M r 
luminosity function based on the method of Blanton 12006). The inset panel 
shows AM,- = - 25 M r — 01 M r , The three filled points indicate the number 
densities of our three z = 0. 1 Main galaxy samples. See Table 1 for the exact 
values of 25 M r for the three 01 M T thresholds. 



estimate is 



CT <7- — 

1=1 



P,' 



P,J 



(2) 



where subscripts i and j denote bins of r p and superscript / 
denotes the subsample. The factor of (N— l)/N corrects for 
the fact that the variance is calculated with samples that are 
slightly smaller than the full sample, thus the variance of the 
subsamples will be larger approximately by the ratio of the 
volumes. The acc uracy of these error e stimates is discussed 
by Z10; see also iNorberg et al.l (120091) . who conclude that 
jackknife errors are generally accurate or conservative on the 
scales investigated here. 

2.2. The maxBCG Cluster Catalog 

The maxBCG algorithm dKoester et al.ll2007bl) utilizes the 
high fraction of early-type galaxies within clusters to identify 
massive halos as overdensities of bright, uniformly red galax- 
ies. The large number of such galaxies within each cluster 
and the tightness of their color distribution allow for accurate 
photometric redshift estimates. Using mock galaxy catalogs, 
the algorithm is optimized to have completeness and purity 
above 90%, or higher as richness increases. The algorithm 
produces a richness estimate for each object, A^oo, that is de- 
fined as the number of red-sequence galaxies brighter than 
25 Mi = -19.2 within an aperture containing a mean overden- 
sity of such galaxies that is 200 times the mean density of 
galaxies. This aperture scales roughly with the v irial radius of 
the ha lo, thus larger halos have larger apertures (lHansen et alJ 
120051) . The richness estimate is not to be confused with the 
occupation number in the HOD; N in the HOD refers to all 
the galaxies (not just red) brighter than a defined magnitude 
threshold within the exact halo radius. The richness A^oOj al- 
though correlated with N to some degree, is used exclusively 
as a reference value with which to bin clusters. Because the 



magnitude limit for maxBCG selection is significantly fainter 
than that of even our M r < -19.5 sample, and the aperture 
within which maxBCG counts member galaxies is not iden- 
tical to the actual halo radius, the values of N are lower than 
those of A200, by roughly a factor of two for 0A M r < -19.5 
and larger factors for brighter thresholds (see Table 2 below; 
also recall that A^oo includes central galaxies while N counts 

only satellites). 

The maxBCG catalog dKoester et alJ l2007ah is a sample 
of clusters identified in 7398 deg 2 of SDSS imaging data, 
roughly the same as the imagin g samples released in DR4 
(lAdelman-McCarthy et al.l [20061) . For our analysis here, we 
utilize a subs e t of th is im aging data that was a nalyzed in 
ISheldon et all d2009bl) and ISheldon et all d2009al) consisting 
of 6325 deg 2 . The redshift range of the clusters in the sam- 
ple is 0.1 < z < 0.3, yielding a nearly volume-limited sam- 
ple of clusters. Although the catalog identifies clusters with 
A^oo > 3, we restrict our analysis to clusters with A^oo > 9. 
Objects with a richness below this value suffer from a higher 
degree of uncertainty due to projection effects. We bin the 
clusters into 6 richness bins listed in Table 2. We note that 
o ur catalog differs slight ly from the publicly released catalog 
in iKoester et alJ d2007al) . which was limited to clusters with 
N200 > 10. 

The average masses of clusters in each richness bin are 
obtained from the w eak gravitational lensing analysis of 
ISheldon et al.1 (l2009bl) . All the clusters in each bin are 
stacked, yielding a high signal-to-noise estimate of the pro- 
jected density contrast profile, AY,(R p ). This profile is de- 
projected and integrated out to a radius at which the mean 
interior density is A = 200/9 c ri t , thereby yielding an estimate 
of the mean cluster mass within each bin. We define this ra- 
dius as ^N200- We differentiate this radius from Z?200c> which 
is the radius for a specific halo, while Rmoo is a quantity spe- 
cific to each richness bin 2 . The weak lensing measurements 
and additional tests with X-ray measurements and velocity 
dispersions of maxBCG clusters demonstrate that A^oo corre- 
lates strongly with dark matt er mass, but there is also a scatter 
between mass and rich ness dBecker et al.ll2007T: iRvkoff et al.l 
l2008t iRozo et ai1l2009l). Th i s scatt er and its uncertainties are 
characterized by IRozo et al.l (120091). 

The cluster masses in ISheldon et al.l d2009bl) are subject to 
a number of biases. First, this analysis assumes that brightest 
cluster galaxy (BCG) identified by the algorithm is located at 
the true bottom of the gravitational potential of the dark matter 
halo. Tests on mock galaxy catalogs reveal that this is not true 
in ~ 10% of observed objects, a fr action that increases m ono- 
tonically with decreasing richness dJohnston et al.l2007h . This 
'miscentering effect' lowers the measured mass relative to the 
true mass. Additionally, errors in the photometric redshifts of 
the sample of background sources behind the clus ters also af- 
fect th e measured halo mass (see the discussion in lRozo et al.l 
120091) . The uncertainties in scatter, miscentering, and weak 
lensing systemat ics are all taken into account in this study. 
In particular, the ISheldon et al.l (l2009bl) we ak lensing masse s 
are estimated to be biased low by 18 ± 6% dRozo et al.ll2009l) . 
We will discuss their quantitative incorporation in the anal- 
ysis in $4] One advantage of our M/N analysis is that it is 
much less sensitive than abundance analysis to miscentering 
errors, since these tend to affect mass and galaxy occupation 

2 Note that i?N 200 is different from the quantity R200. defined in 
lHansen et alJ 120051 ) as the radius at which the mean interior galaxy density 
within clusters is 200 times the mean galaxy density. 



5 



in similar ways. 

2.3. Measuring the Number of Galaxies in Clusters 

To determine the number of satellite galaxies per halo, a 
stacking technique is also used. In each stacked richness 
bin, the total projected galaxy counts are measured in bins 
of log/?,,. To determine the number of galaxies that are as- 
sociated with the cluster and not chance projection, the same 
process of stacking is done with a set of random pointings 
with the same redshift distribution as the clusters. In each 
bin of \ogR p , the mean number of galaxies from the random 
projections is subtracted from the mean number around the 
stacked clusters. Full det ails of this procedur e are given by 
ISheldon et al] (l2009al) and lHansen etaD (120091) . 

Figure QJ shows the projected galaxy number density pro- 
files for the six richness bins in Table 2. These results are 
shown for galaxies brighter than °- 25 M r = -19.0. The density 
profiles do not include the BCG of each cluster, which is lo- 
cated at R p = 0. Thus these data represent satellite galaxies 
only. The error bars are obtained through jackknife resam- 
pling in the plane of the sky with patches of 6.3 deg 2 . To 
obtain the number of galaxies within Rmoo, ^sa\(Rp) is in- 
verted to recover the three-dimensional density profile. The 
inversion is performed using the standard Abel-type integral, 



10 3 - 



(3) 



The projected density is not measured out to infinity, thus in 
practice we fit a power law to the three highest-/?,, data points 
and truncate the integral at 30 h~ l Mpc. Making this upper 
limit twice as large results in negligible differences in the den- 
sity profiles obtained within the cluster radii. Figure QJ) shows 
the inverted three-dimensional density profiles of the satellite 
galaxies. The filled circles indicate the cluster radius in each 
richness bin. These density profiles are integrated out to /?N200 
to determine the total halo occupation of satellite galaxies in 
each richness bin. Figure QJ shows how the number density 
profiles depend on luminosity threshold for a fixed richness 
bin, (A^oo) = 20.8 in this example. As expected, the number 
of galaxies increases monotonically with decreasing luminos- 
ity threshold at all scales. 

The clustering measurements described in § 12. II are for the 
Main sample of SDSS galaxies, which probe z ~ 0.1. The me- 
dian redshift of the maxBCG catalog is z ~ 0.25. For proper 
analysis, we require a consistent sample of galaxies between 
both redshifts. Because the luminosity function evolves be- 
tween these two redshifts, and because 25 M r is not equivalent 
to 01 M r , we choose to define z = 0.25 samples that have the 
same number density as the M ain samples. Figure[2 ]compares 
the l M r luminosity function (Blanton et al. 2003) to our es- 
ti mate of the 25 M r luminosity functi on using the technique 
of lBlantonl d2006l) and emplo yed in the lSheldon et all d2009al) 
analysis. The lBlantonl (120061) technique uses measurements of 
the luminosity function at z ~ 0.1 and z ~ 1 to constrain the 
amount of luminosity evolution and number density evolution 
in both the red sequence and the blue cloud across this redshift 
baseline. These results are used constrain simple star forma- 
tion histories of blue and red galaxies that are used to interpo- 
late the total galaxy luminosity function at z = 0.25. The inset 
panel shows the magnitude shifts between 01 M r and °- 25 M r at 
the three z = 0.1 magnitude thresholds. We will discuss the 
systematic uncertainties of this approach in §2] 
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FIG. 3. — A pedagogical demonstration of the M/N method. Panel (a): The 
circles show the projected correlation function, w p (r p ), measured for o i M r < 
-20.5 galaxies in DR7. The solid curves show HOD model fits to these data; 
in each model, the value of erg is changed. The dotted lines correspond to 
the projected correlation function for the dark matter in each cosmological 
model. Panel (b): The occupation functions inferred from the best-fit models 
to w p (r p ) for each value of og. As <7g decreases, the amplitude of (N)m at 
high masses increases. Panel (c): The points with errors show the maxBCG 
measurements for M/N (now with the miscentering correction from Figure U] 
applied). The curves show the predictions for this quantity for each value of 
trg. For convenience, we have plotted M/N as a function of halo mass rather 
than cluster richness, though in actual model fitting we use richness. 
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FIG. 4. — The effect of cluster miscentering on M/N in mock galaxy cat- 
alogs. The solid curves show the ratio of M/N measured from a sample of 
clusters that are well-matched to halos in the mock catalog to M/N measured 
on the full sample of clusters determined in the mock. Miscentering low- 
ers both the mass and the number of galaxies within clusters, but the effect 
largely cancels in the ratio. The dotted curve connecting the open circles indi- 
cates the effect o f miscentering on the halo masses derived from the AS(i? p ) 
measurements in Johnston et al. ( 2007). In the mock, both M and N are mea- 
sured from the projected mass and number density, respectively, around the 
clusters detected in the mock. We use the mean of all three curves as the mis- 
centering correction on the data, and use the maximal difference of the three 
curves as the error in that correction. 



The error bars on n ga i(r) are correlated, thus to be conserva- 
tive we assume that the errors are 100% correlated; we deter- 
mine the errors on N ssit (< Mj. lm ) by increasing or decreasing 
every « ga i(r) datum by its la error and re-integrating. Our 
"gaiM data are measured on a grid of threshold 25 M r val- 
ues spaced by 0.28 magnitudes. For each threshold, we in- 
tegrate « g ai(/) out to Rmoo to obtain N sat and interpolate to 
obtain N ydt for the exact values of 25 M r that yields the same 
number density as the 0A M r samples (cf. Figure |2}. The val- 
ues and uncertainties of N sdt (< M x / m ) for our three magnitude 
threshold values are listed in Table 2. Figure [TJl shows the 
resulting M/N ratio for our three 0A M r magnitude thresholds. 
Note that the N in M/N here and throughout this work refers 
to the number of satellite galaxies, a choice that makes M/N 
approximately constant from bin to bin, which would not be 
the case if we included central galaxies. To determine the 
statistical error on M/N, we sum in quadrature the fractional 
errors on 7V sa t a nd M f rom the weak lensing measurements in 
ISheldon et ail (l2009bl) . We have 17 total M/N data points, 
6 for the A M r < -19.5 and 0A M r < -20.5 samples, and 5 
for the 01 M r < -21.0 sample. For the brightest sample, the 
richest bin does not have enough clusters to obtain a robust 
density profile. Values of our M and N measurements for the 
maxBCG sample are listed in Table 3. 

When comparing M/N measurements to model predictions, 
we must account for the fact that the M in the M/N ratio 
is the same value for all three magnitude threshold samples. 
Thus the errors are correlated between the samples. We de- 
fine a covariance matrix for the 17 M/N data points where 
i = 1-6 represent the 01 M r < -19.5 sample, ; = 7 - 12 repre- 
sents 0A M r < -20.5, and i =13-17 represents 0A M r < -21. 
Most of the off diagonal elements of the covariance matrix are 



zero, but, for example, data points 1, 7, and 13 are correlated 
because the halo mass is the same for these three data points. 
The diagonal elements of the covariance matrix are obtained 
by total differential error on the M and N data, 



(d{M/N)i\ 2 (d{M/N)i 



(4) 



The off-diagonal terms are also calculated in the same man- 
ner. For example, C\ j is a non-zero element expressed by 



Ci 7 - crl 



d(M/N)i d(M/N) 7 



dM 



dM 



(5) 



where M for the two M/N data points is the same value. We 
note that accounting for this covariance makes little difference 
in our results; the uncertainties in the M/N measurements are 
dominated by systematics that we will discuss in §4. 

As can be seen in Figure [TJl, the M/N ratio is roughly in- 
dependent of cluster richness. In most models of halo oc- 
cupation, the number of satellite galaxies scales as a power- 
law with host halo mass. Most results from observed galaxy 
clustering, as well as N-body simulations, yield a power-law 
index close to unity, which would imply M/N ~ constant 
dKravtsov et al.ll2004t IZheng et alj|2005l; IConrov et al.ll2006t 



van den Bosch et aljT2007b iTinker etall l2007t IZheng et all 
' 20071; lYang et~aT1l2008l) . 



3. THEORETICAL MODELING 

3.1. Galaxy Two-Point Correlation Function 

We parameterize the halo occupation function as two sep- 
arate functions, one representing the occupation of central 
galaxies and one for the occupation of satellite galaxies. For 
central galaxies, we use the standard expression 



cen/M : 



1+erf 



logM-logM^ 

CIobM 



(6) 



where M m - m formally represents the mass at which a halo has a 
50% probability of containing a central galaxy bright enough 
to be contained within the sample. The parameter <7i 0g M is 
related to the scatter in halo mass at fixed luminosity. Func- 
tionally, this parameter controls how "sharp" the transition is 
between halos that host no galaxies and halos that have one 
central galaxy. Because halo mass is monotonically related 
to clustering strength over most of the halo mass spectrum, 
the value of cr\ og M correlates with the large-scale bias of the 
model; a higher scatter brings more low-mass halos into the 
sample, and due to the steepness of the halo mass function 
the resulting large-scale bias is reduced. When M < M*, bias 
is relatively independent of halo mass 3 , thus a\ og M has little 
effect on w p (r p ). For this reason, we allow <7i 0g M to be a free 
parameter for the 01 M r < -20.5 and 01 M r < -21 samples, 
but we fix it at <j\ og M = 0.2 for the faint M r < -19.5 sample. 
In our modeling we adopt a flat prior of 0.05 < o-\ og M < 1.6. 
Values of a\ og M < 0.05 are indistinguishable from cr\ og M = 0, 
while values of a\ og M > 1.6 would be unphysical. 

For satellite galaxies, we adopt an occupation function of 
the form 



(N sat )M = (N cen ) M X 



M 



exp 



-M cl 
M 



(7) 



3 We define the non-linear mass scale M* as the mass at which the linear 
matter variance on the Lagrangian scale of the halo is cr(M) = 1.686. 
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TABLE 1 

SDSS DR7 Spectroscopic Clustering Samples 



0A M r 


o.2S Mr 




Zmnx 


Ng.l 


Hg/Qr 1 Mpc)" 3 


-19.5 


-19.21 


0.068 


0.083 


112497 


1.11 x 10~ 2 


-20.5 


-20.27 


0.104 


0.131 


117588 


3.16 x 10~ 3 


-21.0 


-20.81 


0.126 


0.160 


77381 


1.15 x 10-3 



NOTE. — The number density is corrected for incompleteness. M r is the magnitude 
threshold at z = 0.25 with the same number density as the corresponding 0A M r . 



TABLE 2 
maxBCG Cluster Sample 



N 2 oo range M N2 oo/10 13 h~ l M© N( 0A M r < -19.5) N( 0A M r < -20.5) iV( 0I Afr < -21) miscentering i? N2 oo A" 1 Mpc 



9-11 


3.17±0.12 


4.06 ± 0.08 


1.25 ±0.04 


0.450 ±0.029 


1.33±0.06 


0.59 


12-17 


5.04±0.16 


6.21 ±0.11 


2.00 ±0.06 


0.772 ±0.042 


1.18 ±0.04 


0.70 


18-25 


8.28 ±0.30 


10.40 ±0.22 


3.31 ±0.13 


1.258 ±0.082 


1.14 ±0.04 


0.82 


26-40 


15.54 ±0.59 


17.81 ±0.41 


5.57 ±0.23 


2.26 ±0.16 


1.098 ±0.009 


1.01 


41-70 


24.0± 1.1 


27.92 ±0.85 


8.78 ±0.48 


3.40±0.31 


1.074 ±0.004 


1.17 


71-220 


60.9 ±4.5 


58.0±3.1 


18.8 ±1.7 




1.011 ±0.003 


1.60 



NOTE. — Note: The halo masses are equivalent to those in]Sheldon et al. 1 2009a) with an 18% correction factor as discussed in Rozo et al. ( 2009). 
The N values refer only to satellite galaxies and are given for the equivalent M r thresholds; see text for details. The values listed are the raw 
measurements that do not include the correction for miscentering. The 'miscentering' column indicates the factor, and its error, by which all N 
measurements are multiplied to correct for miscentering. 



Equation (01 parameterizes satellite occupation as a power- 
law at high halo masses with an exponential cutoff at 
low masses, motivated by the results of high-resolution N- 
body and hydrodynamic simulations of galaxy formation 
dKravtsov et al.ll2004t IZheng et alj|2005t IConrovetai] 120061: 
I Wetzel & Whitdl2010l) . The factor of (N cen ) M in equation 
ensures (A^Jm < (N cen )M at all masses. 

We assume that satellite galaxies follow a spatial distribu- 
tion within the dark matter halo of an NFW density profile 
(iNavarro et al.lll997h . However, we do not assume that the 
galaxies trace the dark matter within the halos. The concen- 
tration parameter of NFW profile, c ga i, is a multiple of the con- 
centration parameter of the dark matter, / con = c ga i/c,„. The 
proportionality constant / con is left as a parameter in our anal- 
ysis with a fiat prior of [0.2,2.0], which brackets the extreme 
scenarios for the spatial bias between satellite galaxies and 
dark matter within halos. Thus the shape of the concentration- 
mass relation is the same as the dark matter, but the normal- 
ization is allowed to vary. Fo r the dark ma t ter, w e use the 
concentration -mass relation of iBullock et al.l (1200 lb with the 
parameters of lWechsler et al.l (120061) . 

To calculate the mean number of pairs within each halo, 
we assume that the satellite galaxies are Poisson distributed 
about the mean. This is well supported by both numeri- 
cal re sults (iKravtsov et al.ll2004tlZheng et al.ll2005UGao et al.l 
1201 lb and observational data dLin et al.ll2004l) . Small devia- 
tions from Poisson like t hose shown in the recent resu lts from 
iBusha et alJ (1201 Ol) and iBoylan-Kolchin et alJ ( 1201 Oh do not 
have a significant effect on the clustering in this paper be- 
cause they occur at subhalo masses not probed for the galaxy 
samples analyzed here. 

We use the theoretical model of iTinker et al.l (120051) to cal- 
culate the two-point correlation function in real space, with 
one modification. Because the mass function and bias relation 
used in this analysis are taken from numerical results based on 
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FIG. 5. — The best-fit HOD model fits to w p (r p ) for each 0A M r threshold. 
Uncertainties in the model fits are shown in Figure [6] 



spher ical-overdensity (SO) halo catalogs (ITinker et al.ll2008al 
2010), the halo exclusion must be modified to m atch this halo 
definition. In the SO halo finding algorithm of ITinker et al.l 
halos are allowed to overlap so long as the cen- 



ter of one halo is not contained within the radius of another 
halo. Thus, the minimum separation of two halos with radii 
Ri >R2 is R\, rath er than the sum of the two radii, as done in 
ITinker et al.ld2005l) . 

3.2. Calculating M/N from the HOD 
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Measuring the galaxy content within a cluster sample 
stacked on richness is not identical to measuring the mean 
number of galaxies within the virial radii of halos of mass 
M. First, due to scatter in the mass-richness relation, a stack 
of clusters with the same A200 contains a sample of halos of 
varying mass. Second, each richness bin stacks clusters from 
a range of richness values (e.g., the most massive richness bin 
stacks all clusters with 71 < A^oo < 220). The halo radius 
for each richness bin, Rmoo, is determined such that the mean 
overdensity of all the halos in the bin is 200/9 cl i t . The true radii 
of individual halos in the bin, Rm are not identical due to scat- 
ter; some halos will be smaller the Rmoo while some will be 
larger. Both the scatter and the fixed aperture are accounted 
for when we calculate M/N within our HOD models. 

For each bin of A^oo values, the expected value of M/N for 
a given model is computed by 



M/N'- 



J2 Nl _ m J dMn h (M)P(N 20 o\M)Mf h (R moo \M) 
J dMn h (M)P(N 200 \M)(Ns. dt ) M f g . di (Rmw\M) ' 



18) 



where n i,(M) is the halo mass function given by iTinker et al.l 
d2008al) 4 . P(N2oo\M) is t he probability of a cluster with mass 
M having richness A^oo dRozo et al.ll2009l) . and / is the aper- 
ture correction factor for both the mass and the number of 
galaxies, given the fact that the radius of a halo of mass M, 
Rm, may be different from the fixed aperture /?N200 f° r mat 
richness bin. The integrals in the numerator and denomina- 
tor are evaluated at each value of A^oo in the richness bin and 
summed together. Note however, that /?N200 is me same value 
within a richness bin; see the values listed in Table 2. The 
scatter in richness at fixed halo mass is assumed to be a log- 
normal of the form 



P(N2Qo\M)dN 2 oo = 



1 



exp 



The mean of this distribution is 



2a\ 



A200 



H R = B R +A R ln(M/M } 



pivot J 



(9) 



(10) 



where M pivot = 2.06 x 10 13 x 0. 7 h~ l My, and the f actor of 0.7 
is to change the value listed in iRozo et ail d2009l) from units 
of M Q to h~ l M Q . We will discuss the values and errors on 
the mass-richness scatter in the following section. 

The aperture correction factor is small but non-negligible. 
For a halo with mass M, corresponding to radius Rm, and con- 
centration c(M), the mass enclosed at any radius R is calcu- 
lated by 



f h (R\M) = ^4irp s R 3 y 3 [ln(l + l/y)-(l +yT 1 ] 



(11) 



(c.f., appendix A in iHu & Kra vtsov 2003). where y = r s /R 
and r s = Rm/c{M)). Equation ( fTTT i assumes an NFW form 
for the halo density profile. The parameter p s is a normaliza- 
tion parameter such that fh(RM\M) = 1. Because we assume 
that satellite galaxies also follow an NFW profile, we can use 



4 The ITinker etafl j2008d) mass function is universal at a given redshift 
for constant values of A„,, the overdensity relative to the mean density in the 
universe. The maxBCG observations of mass and Ng^j are made at 200 times 
the critical density (A c ) assuming a value of Q„, = 0.27 and z = 0.25, which 
translates to A m = 508. 



the same scaling relation to determine the number of satellite 
galaxies within an aperture R for a halo of mass M. 

f s dR\M)= — V^^Niooy 3 [infl + l/yMl+yr 1 ] 

{Nsat/M 

(12) 

where p s is once again determined by requiring f gi \(RM\M) = 
1 . As discussed above, the concentrations parameter for satel- 
lite galaxies is defined separately from the dark matter. 

3.3. Probing Cosmology with w p (r p ) and M/N 

To demonstrate the power of combining w p (r p ) and M/N, 
Figure [3] shows both statistics for a series of cosmological 
models. In all panels, the four curves represent four HOD 
models applied to the DR7 clustering measurements of the 
A M r < -20.5 sample, each with different values of a%. For 
demonstration purposes, all other cosmological parameters 
are fixed, and we fix a sa t = 1 ■ For each cosmological model, a 
good fit to the data can be found, demonstrating that cosmol- 
ogy and bias are degenerate for this single statistic. However, 
as the halo population changes with a%, the halo occupation 
must also change in a compensatory fashion in order to obtain 
the same two-point clustering. In Figure[3j:, M/N is shown for 
each of the best-fit HOD models. The low-amplitude cosmol- 
ogy (o"8 = 0.7) requires a larger weighting of high-mass (and 
thus highly-biased) halos in order to match the observed am- 
plitude of the galaxy clustering. This drives M/N lower rela- 
tive to the other cosmologies. The theoretical predictions also 
depend on the value of fi m ; if the matter density increases, 
then the mass of each halo increases proportionately, but the 
number of galaxies per halo remains fixed. Thus O m and og 
will be degenerate in this analysis much as they are in cluster 
abundance constraints 

The comparison here is not one-to-one, given that the data 
are convolved with the mass-richness scatter and that we have 
not marginalized over the many free parameters in our model. 
Nonetheless, this plot demonstrates that the degeneracy of 
cosmology and bias with respect to the two-point correla- 
tion function can be broken by empirical measurements of the 
HOD in the form of M/N. 

4. SYSTEMATIC ERRORS 

Here we detail the many systematic errors that we consider 
in our analysis. Some are associated with theoretical uncer- 
tainties, others with observational uncertainties. 

Uncertainty in the halo mass function: We use the 
simulation-calibrated halo mass function of ITinker et al.l 
d2008al) in all theoretical calculations. This mass function is 
calibrated on spherical overdensity halos, where halo mass 
is defined by the mass within a spherical aperture, the same 
manner as the cluste r mass in th e maxBCG weak-lensing 
mass estimates. The ITinker et al.l mass function also takes 
into account redshift evolution of the mass function, which is 
10-20% bet ween z = and z = .25. Most of the simulations 
analyzed in ITinker et al.l d2008al) are realizations of a single, 
flat-A CDM cosmology re presentative of first-year WMAP re- 
sults dSpergel et al.ll2003h . For this cosmology, the errors on 
the mass function are less than 1% up to nearly 10 15 h~ l M Q . 
For variants around this cosmology, the simulations produce 
a 5% scatter. We therefore incorporate a 5% Gaussian uncer- 
tainty on the amplitude of the halo mass function, introduc- 
ing a marginalization parameter e„ = 1 .0 ± 0.05, which repre- 
sents the amplitude of the halo mass function relative to the 
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FIG. 6. — Each panel shows the fraction difference between the measure- 
ments of w p (r p ) and the best-fit model for each o i M r threshold. The contours 
indicate the range of models from the MCMC chain. The inner contours in- 
clude 68% of all elements in the chain, while the outer contours bracket 95% 
of the elements in the chain. The data points are highly correlated, thus we 
caution the reader against x-by-eye. Note that all models are simultaneously 
fit to M/N as well as w p (r p ). 



iTinker et al.l d2008al) fit. We describe the incorporation of this 
and subsequent uncertainties at the end of this section. 

Uncertain ty in the large-scale halo bias relation: The 
mass f unction is coupled to the bias functions of 
d2010l) . These bias functions are calibrated on the 



Tinker et al 



Tinker et al 



same N-body simulations as the mass function and with the 
same halo definitions. Thus, the abundance and clustering of 
halo s are always calculat ed using a self-consistent halo defini- 
tion. ITinker et al.l (1201 Oh find a 6% scatter in bias among the 
simulations. To be conservative, we implement a 6% Gaus- 
sian error on the amplitude of the halo bias relation. As with 
the mass function, we delineate the prior on the bias ampli- 
tude as Eb, with a mean of 1.0 and an uncertainty 0.06. 

Uncertainty in the scale-dependence of the halo bias: At 
scales of r < 5 h~ l Mpc, halo clustering deviates from a scale- 
independent factor of the matter clustering. This scale de- 
pendence is minimized when the non-linear matter clustering 
(rather than the linear) is used to define the halo bias. We 
as sume that the shap e of the scale-dependence follows that 
of ITinker et al.l (120051) but allow the magnitude to differ. We 
include a 15% Gaussian error in the deviation from constant 
bias, i.e., 

Sb = b(r)-b , (13) 

where b(r) is the scale-dependent bias of ITinker et al.l (120051) 
(with a modification for very small scale halo clustering 
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FIG. 7. — Each panel shows the M/N measurements, normalized by 
Pent/fig, for each threshold sample. The solid curve in each panel shows 
the best-fit HOD+cosmological model. The inner and outer shaded regions 
bracket the 68% and 95% ranges of models in the MCMC chain, respectively. 
Note that all models are simultaneously fit to w p (r p ) as well as M/N. 



given in Appendi x A) and bo is t he large-scale bias of 
ITinker et al.ld2010l) . The ITinker etaTI (120051) scale-dependent 
bias, b(r)/bo, asymptotes to unity at r > 5 /z _1 Mpc. Param- 
eterizing the uncertainty in the scale dependence through eq. 
( fT3] l ensures that bias becomes linear at large scales. We refer 
to this prior as e,-, with an uncertainty of 0. 15. 

Uncertainty in the par ameters of the mas s-observable rela- 
tion: As described above. IRozo et al.l(l2009l) describe the rela- 
tion between halo mass and observed cluster richness from the 
maxBCG algorithm as a power law with parameters o.r and 
f3 R . The scatter of mass at fixed richness is modeled as a log- 
normal with var iance erg. Using X -ray data from the ROSAT 
All-Sky Survey. IRozo et al.l (1201 Ol) both calibrate this relation 
and determine the uncertainties in each parameter. These are 
found to be 



A R = 0.750 ±0.024 
B R = -1.09 ±0.09 
cr R = 0.35 ±0.07. 

These three parameters are also added to the chain with the 
uncertainties listed above. 

Uncertainty in the evolution of the HOD: Due to the redshift 
baseline of our two galaxy samples, z = 0.1 and z = 0.25, we 
must account for possi ble evolution in the HOD over this nar- 
row range. Analysis by lZheng et al.l(l2007l) over a much larger 
baseline, from z = 0.1 to z = 1 using the clustering of DEEP2 
galaxies, demonstrates that the M sat /M m i n ratio evolves from 
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FIG. 8. — Constraints in the Q,„-ag plane, marginalizing over all other 
parameters and applying the priors listed in Table 3. The blue contours show 
68% and 95% constraints from the M/N results. The yellow and grey shaded 
region indicates the lcr and 2a constraints, respectively, from ITinker et al.l 
120051) . The green contour sho ws the constraints from WMAP7 (CMB alone, 
assuming a fiat-ACDM model;. Komatsu et al. 20 1 (J) . The red contours show 
the combined constraints from M/N and WMAP7. 



18 to 16 (independent of luminosity). lAbbas et all (1201 Oh 
see even less evolution from z = 0.4 to z = 1.1, but using the 
smaller zCOSMOS data set. 

We have investigated this issue theoretically in two ways: 
through analysis of the Mill ennium semi-analytic galaxy cat- 
alogs of iBower et al.l (120061) , and through the hi gh-resolution 
N-body simulation Bolshoi dKlypin et alJl2010h . Recall that 
we are analyzing samples constructed to have the same space 
densities at their respective redshifts, not the same magni- 
tudes. This minimizes the effect of the evolution of the galaxy 
population. In the semi-analytic Millennium model, the HOD 
at fixed fig shows a moderate increase of 5-10% in (N^m at 
a given M from z = 0.1 to z = 0.25. 

In the Bolshoi simulation, we investigate the the evolu- 
tion of the halo occupation of subhalos. The connection be- 
tween substructure and galaxies is well-established (see, e.g. , 



Kravtsov et al 



Behroozi et al 



2001 IConrov"eTaTi r2006: Mo ster et al.ll20ia 



20101 for some examples), thus this method 



can be used to construct galaxy catalogs without attempting to 
model galaxy formation physics. At fixed space density 5 , the 
halo occupation of subhalos shows < 10% evolution over our 
redshift range (Reddick et a l, in preparation), va rying slightly 
with galaxy number density. ISimha et~aT1 d2010t) find that sub- 
halo occupation statistics track the the statistics of galaxy oc- 
cupations in smoothed particle hydrodynamics (SPH) simula- 
tions fairly accurately, with the largest deviations arising be- 
cause a small but not negligible fraction of SPH galaxies suf- 
fer dramatic stellar mass loss after entering high mass halos, 
thus ending up less massive than predicted. 

5 Here we mean density of all halos and subhalos above a given maximum 
circular velocity. For subhalos, the value of V max at the time of accretion is 
used to connect to galaxy luminosity. 



Based on these results, we incorporate HOD evolution un- 
certainty as a 10% Gaussian multiplicative error on M/N, 
centered on no evolution. 

Uncertainty in the luminosity function at z = 0.25: As odd 
as it may seem, there is no reliable measurement of the 25 M r 
luminosity function. As discussed above and shown in Figur e 
12 we infer this function by the method of Blantonl d2006l) . 
Although the redshift difference between the Main sample 
and the maxBCG sample is relatively small, the amplitude 
of M/N is sensitive to the details of our luminosity function 
estimate. The key quantities are the 25 M r magnitudes that 
yield the same number densities as the z = 0.1 samples. A 0.1 
mag error in these values yields a ~ 20% error in M/N. Given 
that the overall evolution of the luminosity function is 0.2 to 
0.3 magnitudes, we use 0. 1 magnitudes as a conservative error 
estimate on our estimate of the z = 0.25 luminosity function. 
As with the evolution of the HOD, evolution in the luminos- 
ity function shifts the amplitude of M/N, thus we include an 
additional 20% Gaussian multiplicative error on our theoret- 
ical calculation of M/N for a given model. The quantitative 
incorporation of this uncertainty with HOD evolution will be 
discussed at the end of this section. 

Uncertainty in the calibration of the weak tensing masses: 
iRozo et al.1 d2009l) apply a correction factor of 1 .18 ± 0.06 to 
the maxBCG halo ma sses inferred through wea k lensing in 
ISheldon et all (l2009bl) and lJohnston et al.l J2007I) . This factor 
is to compensate for the fraction of sources with redshift er- 
rors large enough to place their true redshift in front of the 
lensing clusters, thereby diluting the shear signal. An error in 
halo mass scale is partly compensated by an increase in the ra- 
dius at which the mean interior density is 200/9 cl -j t , which is the 
radius at which N sat is measured. When measuring M/N for 
both the original ISheldon et al.1 (l2009bl) masses and the new 
corrected masses, we find that the 18% increase in the mass 
at each richness bin yields only a 9% increase in the M/N ra- 
tio at each richness bin. The remaining 6% uncertainty in the 
halo masses therefore yields a 3% uncertainty in M/N. This 
error is correlated over all data points, and it is included with 
the other data systematics listed above: evolution of the HOD 
and the luminosity function. 

Miscentering of the maxBCG algorithm: As with all opti- 
cal cluster finders, maxBCG is not 100% pure and complete. 
Some fraction of the time, the correct halo is identified but the 
wrong galaxy is identified as the BCG Additionally, chance 
projection of two clusters can result in the two objects being 
m erged in the resulting c luster catalog. In the mass estimates 
of ISheldon et alJ (l2009bl) . these effects are not accounted for. 
They are also not accounted for in the measurements of the 
galaxy counts. Because satellite galaxies follow roughly the 
same distribution as the dark matter within the cluster, the ef- 
fect on their ratio is significantly sm aller than on either quan- 
tity individually. IRozo et all d2010h find that the fraction of 
chance projections is small compared to the intrinsic scatter 
in the mass-richness relation, thus they conclude that the ef- 
fect of projections on their constraints is negligible. 

We quantify the effect of miscentering using mock galaxy 
catalogs on which the maxBCG algorithm is run. The mock 
catalogs use the ADDGALS (Adding Density Determined 
Galaxies to Lightcone Simulations) algorithm to pop ulate a 
dark matter simulation with galaxies dWechsl er 2004; Wech- 
sler et al preparation). Galaxies are assigned to dark matter 
particles based on their dark matter densities, constrained to 
match the observed two-point galaxy clustering statistics and 
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luminosity function, and are assigned galaxy colors based on 
the observed distribution of galaxy color at a given galaxy 
density and luminosity. Mocks created with this algorithm 
has been used to test the maxBCG algorithm in several pre- 
vious works (e.g . iKoester et al.ll2007at I Johnston et al.ll2007t 
iRozo et al]l20ia iGerdes et al.ll2010l) . Clusters in which the 
detected BCG is the same as a central galaxy in a mock dark 
matter halo are deemed "well-matched". Clusters where this 
match could not be made represent objects where the detected 
BCG is not a central galaxy. For the well-matched sample we 
repeat all the steps described in Figure Q] For each bin in 
richness, the mock background-subtracted surface density of 
galaxies is measured and inverted to obtain n ga \(r). This quan- 
tity is integrated out to /?haio, yielding N ydt . The same process 
is performed on the dark matter particles as well to mimic 
the effect of miscentering on the weak lensing masses and the 



M/N ratio. To obtain Mh a i and /?haio the dark matter den- 
sity profile is integrated out until A = 200/9 cr j t . The results are 
shown in Figure|4]for all three luminosity thresholds. At high 
richness, miscentering has a negligible effect on M/N. At 
Afcoo ;$ 20, M/N for the well-matched sample is larger than 
that of the full sample of clusters. This implies that the effect 
of miscentering on mass is somewhat stronger than on galaxy 
number. However, the effect of miscentering on M/N is sig- 
nificantly smaller than on M by itself, as shown by the circles 
connected by the dotted line. Although there appears to be 
a monotonic effect with luminosity at low richness, it is not 
clear if this trend is real or noise. We take the mean of these 
three curves as the miscentering correction, with the range 
between the curves as a ler Gaussian error on the correction. 
This error is added in quadrature to the statistical errors on 
M/N. It is listed in Table 2. 
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Environmental dependence of the HOD: It has been demon- 
strated that the properties of halos at fixed mass depend on 
large-scal e environment, an effect dubbed 'assemb ly bias' 
(see, e.g.. ETao & White! 120061: IWechsler et al.ll2006l) . Semi- 
analytic models of galaxy formation have predicted that this 
assembly bias could propagate into the properties of galax- 
ies within halos at fixed mass (ICroton et al.ll2007t IZhu et all 



120061) . although results from hydrodynamic simulat ions do 
not support this (iBerlind et al.ll2003l; lYoo et al.ll2006b . In the 
standard HOD, (N)m depends only on M and not on a second 
parameter, such as 8. If occupation depends on environment, 
then the large-scale bias of a galaxy sample may differ from 
that calculated through the standard HOD, biasing the results 
of the modeling. By modeling both w p (r p ) and the distribu- 
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tion of galaxy voids in the SDSS. lTinker et al.l d2008bl) found 
that any assembly bias must have less than 5% effect on the 
large-scale clustering of luminosity-defined samples for the 
luminosity range probed by SDSS data. Furthermore, our ex- 
clusion of clustering data at r p > 3 h Mpc attenuates any im- 
pact assembly bias may have on our HOD fits. We conclude, 
given the small effect of assembly bias and the exclusion of 
large-scale data, that our 10% prior on the evolution of the 
HOD subsumes any error accrued due to assembly bias. 

Combination: The uncertainties in evolution of the HOD, 
the z = 0.25 luminosity function, and the calibration of the 
weak lensing masses all have the effect of systematically 
shifting the M/N measurements up or down proportionately 
for all richness bins. These three systematic errors are uncor- 
rected, thus we include a \/l0 2 + 20 2 + 3 2 = 22.6% systematic 
error on all the M/N data. To incorporate this uncertainty 
in the MCMC analysis, we add a 'systematic bias parameter' 
/ sys to the chain. For each element in the chain, the HOD pre- 



diction for M/N is multiplied by / sys . We enforce a Gaussian 
prior of 1 .0 ± 0.226 on / sys . 

5. RESULTS 
5.1. Fitting the Data 

We determine constraints on all free parameters through the 
Monte Carlo Markov Chain (MCMC) method. For each ele- 
ment in the chain, we determine the x 2 for the model in the 
following manner. First, x 2 is calculated independently for 
all three sets of w p (r p ) measurements using the full covari- 
ance matrix for each sample. Next, the x 2 for the 17 M/N 
data points is calculated using eq. (0 and the covariance ma- 
trix. The total x 2 for each model is the sum of the two. To 
marginalize over the systematic uncertainty in the mass func- 
tion, the parameter e„ is a free parameter within the chain. For 
each element i in the chain, the mass function for that model 
is nhj{M) = n/,(M) x e„. While running the chain, we adopt a 
5% Gaussian prior on this parameter. We marginalize over all 
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TABLE 3 
List of Parameters 



Name 


Type 


Value 


Meaning 


1 Jlyf 


observational 




SDSS r-band absolute magnitude k+e corrected to z — 0.1. 




observational 


— 


SDSS r-band absolute magnitude k+e corrected to z = 0.25. 




cosmological 


no prior 


matter density relative to the critical density 


CT8 


cosmological 


no prior 


amplitude of linear matter fluctuations on the 8 h Mpc scale 


ho 


cosmological 


0.72 ±0.05 


Hubble constant in units of lOOkm/s/Mpc 


n s 


cosmological 


0.96 ±0.03 


spectral index of primordial fluctuations 


^min 


HOD 


no prior 


central occupation function; see eq. ((6). 


°"logM 


HOD 


[0.05,1.6] 


halo mass-galaxy luminosity scatter; see eq. (6). 


M sM 


HOD 


no prior 


amplitude of the satellite occupation function; see eq. Q. 


M cut 


HOD 


no prior 


cutoff of the satellite occupation function; see eq. (7). 


Osat 


HOD 


no prior 


power-law slope of the satellite occupation function; see eq. |7). 


/con 


HOD 


[0.2,2] 


satellite galaxy density profile: / con = Cgal/ c halo 


e n 


systematic 


1.00 ±0.05 


uncertainty in amplitude of halo mass function 


c-b 


systematic 


1.00 ±0.06 


uncertainty in amplitude of halo bias relation 


e r 


systematic 


1.00 ±0.15 


uncertainty in amplitude of scale-dependent bias 


/sys 


systematic 


1.00 ±0.22 


uncertainty in the measured amplitude of M/N 


Br 


maxBCG 


-1.09 ±0.09 


amplitude of mass-richness relation; see ecr jlOl 


Ar 


maxBCG 


0.750 ±0.024 


slope of mass-richness relation; see eq. ill 01 


vr 


maxBCG 


0.35 ±0.07 


scatter in mass-richness relation; see eq. l|9} 


N200 


maxBCG (observed) 




richness of maxBCG clusters 




maxBCG (observed) 




cluster radius at which A = 200p cr jt in a given richness bin 


W S at 


maxBCG (observed) 




true number of satellite galaxies in maxBCG clusters 



NOTE. — The HOD parameters, M m i n , M sa t, «sai, Mcut, /con, are constrained separately for all three clustering samples. Thus 
there are three values of each parameter that included in the analysis. The HOD parameters <Ti ogM is free for the brighter two samples 
and fixed at <ti os m = 0.2 for the faint sample. M mm for each sample is set by the number density of each sample once the other HOD 
parameters are specified. Thus there are 14 free HOD parameters total in this analysis. 




other sources of systematic uncertainty in the same manner, 
i.e., we repeat this procedure for e/,,e,-,/ sys , and the parame- 
ters of the mass-richness relation. 

Our model has a total of 21 parameters and priors. A list 
can be found in Table 3. There are 14 HOD parameters; each 
galaxy sample is modeled with an HOD of 4 or 5 parame- 
ters (we exclude a\ og M from the l M r < -19.5 sample). There 
are two free cosmological parameters, fl m and a%. There are 
two additional cosmological priors, /io and «„ which enter 
into the calculation of the linear matter power spectrum for 



a given cos mological model , for w hich we use the transfer 
function of lEisenstein & Hul d!999t) . We enforce Gaussian 
priors on these parameters that are somewhat broader than 
their current uncertain ties considering the latest CMB results 
dKomatsu et alJl20Toh . As we will demonstrate, the influence 
of these parameters is minimal, but we include them for the 
purpose of marginalizing over acceptable shapes of the matter 
power spectrum. There are 7 additional priors on our system- 
atic errors; 3 for dark matter halo statistics, 3 for the clus- 
ter mass-richness relation, and one — / sys — to incorporate sys- 
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tematic uncertainties in our measurements. 

The x 2 for the best-fit model is 26.9. With 16 free param- 
eters for 38 data points, this model yields a x 2 P er degree of 
freedom of 1.22. The probability of drawing a value > 26.9 
from a x 2 distribution with 22 degrees of freedom is 0.215. 
Figure [5] shows the projected correlation function measured 
from DR7 and the best-fit model. Due to the vanishingly small 
error bars, we also show the fractional differences between the 
model and data in Figure [6] The shaded bands show the 68% 
and 95% ranges in the models from the chain at each value of 
r p . Figure [7] presents the maxBCG M/N measurements and 
the results from the MCMC chains. In each panel, the curve 
shows the best-fit model, while the shaded regions once again 
indicate the 68% and 95% ranges at each value of A^oo- The 
breakdown of the total x 2 by the different data sets is as fol- 
lows: x 2 =61, 6.5, and 5.7 for the l M r < -19.5, -20.5, and 
-21.0 clustering samples, respectively. For the M/N data the 
X 2 = 10.6. Parameter values and uncertainties for the HOD pa- 
rameters are given in Table 4, while cosmological constraints 
under various priors and assumptions are listed in Table 5. We 
note that the best-fit value of / sys is 1.13 — given our model 
assumptions, fitting the data favors M/N values 18% higher 
than our direct estimates, but within our estimated systematic 
uncertainty. 

5.2. Cosmological Parameter Constraints 

Figure [8] shows the main cosmological constraints in the 
plane. In cluster abundance studies, there is a natural 
"banana curve" degeneracy between fl,„ and erg; the number 
of clusters above a fixed mass varies with the overall matter 
density, while the number of massive objects is highly sensi- 
tive to the amplitude of matter fluctuations. This degeneracy 
curve also exists for cluster M/L or M/N ratios, for the rea- 
sons described in $3.31 T he shaded banana c urve in Figure 
[S] shows the results of the iTinker et al.l (120051) analysis. The 
degeneracy between Cl m and erg is still seen in our analysis, 
with a combined constraint of fi„ = 0.465 ± 0.026. The 
power on f2 m was chosen such that the likelihood of the pa- 
rameter combination was symmetric and yielded the small- 
est fractional error. Without adding any additional datasets 



and marginalizing over all other parameters, the constraints 
on individual cosmological parameters are: f2,„ = 0.29 ± 0.03, 
rj8 = 0.85 ±0.06 (68%). We emphasize that these constraints 
follow from SDSS galaxy clustering and maxBCG galaxy 
and weak lensing profiles, given only loose priors on the pri- 
mordial power spectrum and the broad assumptions of our 
HOD framework. Combining these results with the latest 
constraints f rom WMAP7 (CMB- only constraints, assuming 
flat ACDM; iKomatsu et alJ[20loh . our constraints are Q m = 
0.280 ±0.012, e rg = 0.812 ± 016. 

Although the ITinker et al.l (120051) result overlaps with our 
ler error contour in the Sl„,-erg plane, the best-fit value of the 
cluster normalization is higher in our current results. Thi s off- 
set is driven primarily by the differences in the amplitude of 
the M/N measurements from the maxBCG catalog relative to 
the mean M/L ratio of the 17 clusters in the CNOC2 survey. 
However, there are several marked differences in the present 
analysis. First, we model M/N as a function of cluster rich- 
ness, as opposed t o using the mean M/L ratio. Second, in 
ITinker et al.l (120051) we fixed a sat = 1 . In this paper, we find 
1 .0 < a sa t < 1 • 1 depending on luminosity. If we had fixed 
a sat = 11 in the previous analysis, this would yield a 6% 
increa se in the ^ m ' 6 erg constraint (see Fig. 4 in ITinker et al.1 
120051) . Lastly, systematic errors in the HOD modeling, as well 
as redshift evolution fro m the mean CNOC2 redshift z ~ 0.3, 
were not included in the ITinker et al.l (120051) analysis. This 
would increase the errors and bring the two analyses into bet- 
ter statistical agreement. 

Figure|9]illustrates the origin of the constraints in the £l m -cr% 
plan, using data from the 0A M r < -20.5 sample. We choose 
four models from the MCMC chain that lie on our 2cr contour, 
two that bracket the width of the Sl m ' 5 erg degeneracy axis (both 
erg and f2 m are low or high) and two that bracket the range of 
allowed models along a constant il[J; 5 erg. Figure [2) plots that 
fractional discrepancies of each model with the w p (r p ) data, 
while Figure [5J: compares the predicted and observed M/N. 
The M/N agreement appears acceptable in all cases, but that 
is because the M/N normalization systematic / sys is allowed 
to vary. For the high and low Sl^Vg models, / sys has been 
forced to unlikely values that account for much of the % 2 . Fig- 
ure 9d shows the M/N predictions when / sys is fixed to our 
best-fit value of 1.13. In this case, the high and low f^'Vg 
models predict M/N values that are clearly too high and too 
low, respectively. The two models that define the allowed 
range along the degeneracy axes are ruled out largely by their 
discrepancies with the clustering data, especially for the low- 
rim, high-erg model. Somewhat counter-intuitively, the dis- 
crepancies with w p (r p ) have the same sign in both cases, but 
the M/N data are also constraining the HOD — in the low-erg 
model, lowering (A^Jm at high masses would reduce w p (r p ), 
but it would also raise M/N values are already too high. 

Figure [10] demonstrates the effect of systematic uncertain- 
ties on our cosmological constraints. In the first panel, our 
uncertainties on the halo mass function, halo bias relation, 
and scale dependent bias, are reduced to 1%, 1% and 3%, 
respectively, which dramatically shrinks the allowed param- 
eter range along the £l m ' 5 crg degeneracy axis. The tightening 
of the constraints is due, in order of importance, to the re- 
duced uncertainty in the large-scale bias, the scale-dependent 
bias, and the mass function. The uncertainty in the mass func- 
tion has nearly negligible effect on the results, which demon- 
strates once again the distinction between M/N and cluster 
abundance constraints. The middle panel shows the effect of 
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the uncertainty on the systematics in the measurements — the 
evolution of the luminosity function, the HOD, and the weak 
lensing calibration, encapsulated by / sys . Here we reduce the 
combined uncertainties to 1%, which leaves the range along 
the degeneracy curve unchanged but shrinks the width of the 
£l m -a& degeneracy curve by about a factor of two. We note 
that once the prior on / sys goes below 5%, the results are un- 
changed. At that level of precision, we are statistically limited 
by the M/N signal measured from the cluster sample. The far 
right panel shows the cosmological constraints when reduc- 
ing both the theoretical and measurement uncertainties. For 
a given galaxy clustering measurement, the M/N data set the 
amplitude of the m -(7g degeneracy axis. The value of / sys 
varies the amplitude of M/N, thus widening the error con- 
tours parallel to this axis. The shape and amplitude of the 
galaxy clustering measurements constrain the length of the 
error contours along the axis. By reducing the systematics, 
it would be possible to tighten the constraints on a% and VL m 
to ~ 3-5% from this technique alone, without any additional 
data sets (see the lower three rows in Table 5). 

Figure QT| explores the effect of other priors in our analy- 
sis. In the left panel, we enforce a 10% Gaussian prior on 
the mass-concentration relation for satellite galaxies. This 
has minimal impact on the cosmological constraints, demon- 
strating the insensitivity of this approach to the details of the 
spatial bias of galaxies within clusters — our method is only 
sensitive to the mean number of satellites. The middle panel 
shows the effect on tt m -(7& if we tighten our priors on ho and 
n s , which has the effect of narrowing the range of shapes of 
the linear matter power spectrum. Once again, there is little 
effect on cosmology from the shape of P{k) — our method is 
most sensitive to its amplitude. The right panel assumes min- 
imal error on the mass-richness relation and its scatter. Here 
again the change in the cosmological constraints is minimal. 

Figure Q~2] shows the constraints on cosmological parame- 
ters as well as the marginalization parameters. In these con- 
tour plots one can see clearly the dependence of our cosmo- 
logical parameters on / sys , q,, and e r , as well as the lack of 
correlation with e n , ho, and n s . Of our systematic uncertain- 
ties, Eb (the normalization of the halo bias relation) shows 
the strongest degeneracies with the cosmological parameters 
f2 m and og- The uncertainty on halo bias anticorrelates with 
<Tjj. The clustering data determines the parameter combina- 
tion b sa i<Ts, thus the same fit to w p (r p ) can be achieved with 
a lower erg and an artificially enhanced bias factor. Because 
the M/N data determine the Cl^fa% combination, ej shows a 
positive correlation with O m to preserve the cluster normal- 
ization as erg varies. The red curves plotted on top of the his- 
tograms represent the prior applied to that parameter. The his- 
tograms for ho, n s , eb, and e n are coincident, or nearly so, with 
their priors. The best-fit marginalized values of e r and / sys are 
~ -1 .3(7 and +0.6(7 off their prior values, respectively. For e, , 
a value below unity brings the scale-dependent bias formula 
closer to scale-independent. At scales below ~ 1 h~ l Mpc, 
this increases the amplitude of the halo-halo clustering. The 
value of e r is positively correlated with ag through their effect 
on w p (r p ) at the transition scale between pairs of galaxies in 
distinct halos and pairs of galaxies that exist within a single 
halo (commonly referred to as the 'one-two to two-halo tran- 
sition region'). Pairs of galaxies at the transition regime can 
come from both sources — satellite galaxies in massive halos 
and galaxy-galaxy pairs in smaller halos (that can be closer 
together). If two-halo clustering is higher at r p < 1 h~ l Mpc, 
there is less need for massive clusters to supply the pairs to fit 



w p (r p ) at that scale. The offset between the MCMC results for 
/ sys and its prior indicate some tension between the clustering 
measurements and the maxBCG M/N data. Recall that / sys 
multiplies the model calculations of M/N; to fit the data with 
/ sys = 1 requires a higher a% than is preferred by the clustering 
measurements. However this offset is within our uncertainty 
of the evolution between the redshifts of the two samples. Al- 
though circumstantial, it is worth noting that the best fit value 
of /sys is consistent with t he amount and direc tion of evolu- 
tion in the HOD seen in the lBower et alj d2006l) semi-analytic 
galaxy formation model (see Appendix B), and with the di- 
rection of the HOD evolution seen in our abundance matching 
models. 

5.3. Halo Occupation Constraints 

The lcr parameter constraints for all three clustering sam- 
ples are listed in Table 4. Figure Qj] shows the constraints on 
the HOD parameters for the A M r < -19.5 clustering sam- 
ple. There is a strong degeneracy between a sat and M sat . As 
M sat increases, the number of galaxies in high mass halos de- 
creases, but this can be compensated for by increasing a sat . 
Tilting the (iV S at) m power law has the effect of reducing the 
number of satellite galaxies in halos with M < M sat . It is satel- 
lite galaxies in these low-mass halos that provide the pairs in 
the correlation function at the smallest scales probed by our 
w p (r p ) measurement, r p < 0.3 /j _1 Mpc. Thus to fit w p (r p ), 
when a sat and M sat increase, M cut decreases. In the best-fit 
model for this sample, satellite galaxies are more centrally 
concentrated than the dark matter by a factor of 1.6, but the 
constraints on this parameter are weak and not strongly corre- 
lated with any other HOD parameter. 

For 01 M r < -20.5, shown in Figure [14] similar degenera- 
cies are seen between a sat , M sat , and M cut . There is a tail in the 
likelihood function for M cut to low masses; when M cut < M m ; n , 
M cut no longer has an effect on the HOD or w p (r p ), thus all 
values below this mass scale are equally likely. Rather than 
place a prior on M cut that depends on the value of M m ; n , we 
allow the chain to roam free through this part of parameter 
space, which yields unusual likelihood shapes. For this sam- 
ple, values of a\ og M > 1 are strongly excluded because larger 
values of the scatter decrease the large-scale bias of the sam- 
ple and cannot be reconciled with the observations even given 
the range of erg values probed in the MCMC chain. For this 
sample, the spatial distribution of the satellite galaxies is best 
represented by the concentration-mass relation adopted for 
the dark matter halos, but the constraints on / con are also weak 
and do not depend on the other HOD parameters. 

For the bright sample, shown in Figure [14] cr\ og M is more 
tightly constrained; low values of the scatter produce too high 
a large-scale bias relative to the w p (r p ) measurements. A 
value of cr\ oM that monotonically increases with galaxy lumi- 
nosity is expected if the scatter in galaxy luminosity affixed 
halo mass is a constant, as widely assumed. The strange like- 
lihood function for M cut is due to the fact that the M cut < M m j n 
for most of the chain, as discussed above. The / con parameter 
is unconstrained in this model; the error bars in the one-halo 
term are somewhat larger because of the smaller galaxy num- 
ber density and the lower fraction of galaxies that are satellites 
relative to fainter samples (see, e.g., Z10). 

For all three samples, the power-law index on the satellite 
occupation function, a 8a t, is close to unity, as expected from 
theoretical studies and from previous analyses of SDSS clus- 
tering. There is some tension between our values of a sa t and 
those derived from the clustering-only analysis of Z10. For 
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TABLE 4 
HOD Parameters Constraints 



Name 


01 M r < -19.5 


0A M r < -20.5 


01 M r < -21.5 


logM min 


11.59 ±0.07 


12.21 ±0.11 


12.87±0.12 


logM sat 


12.94 ±0.06 


13.46 ±0.05 


13.87 ±0.05 


a sa t 


1.01 ±0.04 


1.03 ±0.05 


1.08 ±0.05 


logM cllt 

log<7l ogM 
/con 
/sal 


12.48^;!! 
-0.69 
1 61 +ft34 

1 - D1 -0.44 

0.16±0.02 


12 60 +0 25 

_n =4+0.40 

n o 7 +0.45 
u - y '-0.23 

0.15±0.01 


1 29+ - 38 

-o.i2^;8l 

0.81^ 
0.13±0.01 



NOTE. — All halo masses are in units of hr l Mq. All logarithms are in base-10. / sa t is the 
fraction of galaxies that are satellites. We note that any comparison between these results and 
others should account for any difference in halo mass definition. 



TABLE 5 

Cosmological Parameter Constraints 



Name Q„, <r s ^S; 5 °"8 

M/N 0.29 ±0.03 0.85 ±0.06 0.465 ±0.026 

M/A'+counts 0.280^020 0.809^^ 

M//V+WMAP7 0.290 ±0.016 0.826 ±6.020 

M/JV+counts+WMAP7 0.280 ±0.012 0.812 ±0.016 

M//V+WMAP7+BAO+// 0.282 ±0.011 0.8 19 ±0.015 



C8 

M/N minus theory sys. 0.276^ ^ 0.908_^ g^ 

M/N minus data sys. 0.331;*'' 2 !? 0.842^ : ^o 0.493^^ 

M/N minus all sys. 0.297;™!]' 0.952™ : <f * 

NOTE. — Results including "counts" refer to the res ults fr om Rozo et al. 1 201,3). Results 
using WMAP7 data refer to seven-year CMB data from Komatsu e t al.l 120 lOl) . These results 
assume a flat ACDM cosmology. WMAP7 and WMAP7+B A0+f/p results are t aken from the 
publicly available MCMC chains from http : / /lamb da ■ gsf c ■ nasa . gov/ The bottom 
thre e parameter sets show constraints by reducing the errors on our systematics, as in Figure 
[lO](no additional data are included). 'Minus theory' means setting the priors on e„, e/,, and e r 
to 1%, H 

all' means combining both theory and data results. 



01 M r < -20.5, and 0A M r < -21 .0, Z10 find a sat = 1 .15 ±0.03 
and a S at = 1.15 ±0.06, respectively. Our analysis yields 
«sat = l-03±g$ and a sat = 1.08±0.05. Some of this difference 
may be due to the choice of erg = 0.8 in Z10; this is la be- 
low our best-fit value. But to account for the full difference, 
the M/N data must be included. It is notable that our con- 
straints on HOD parameters are not appreciably better than 
those of Z10 even with the addition of extra data on the HOD; 
however, Z10 fix their cosmology, which naturally tightens 
the constraints on the HOD parameters, and the include larger 
scale w p (r p ) data. Z10 use a cutoff power-law to parameterize 
(N s1 h)m, as opposed to the exponential cutoff in eq. (|7), and 
they assume Ahak> = 200p„, for all calculations, thus a one-to- 
one comparison of the two is not straightforward. 

6. DISCUSSION 

6.1. Comparison to Cluster Abundance Studies 

Throughout the analysis in this paper, the information from 
cluster number counts — i.e., the space density of clusters as 
a function of richness or mass — has not been used. Figure 
[l6]compares our results to those of various cluster abundance 
studies. In particular, the left-hand side compares the cos- 
mological constrai nts from M/N to th ose of abundances from 
maxBCG clusters dRozo et al.ll2010l) . Both of these analy- 



ses use the same sample of clusters, but in nearly comple- 
mentary ways. The M/N results are both in good agreement 
with the abundance constraints and competitive with the abun- 
dance constraints. The right-hand side of Figure [TSIpresents 
the 95% confidence region in the fl m -o-& plane for the two 
maxBCG studies in comparison to several cluster abundance 
results from X-ray observations. Here again, the M/N results 
are in good agreement with and competitive with these ap- 
proaches. 

There are several ways in which the M/N approach is 
complementary to and independent of the cluster abundance 
method, even when using the same cluster sample. The pa- 
rameters e r , et,, and the evolution of the luminosity function 
have no impact on abundance analysis. The dominant the- 
oretical systematic uncertainty for cluster abundances is the 
shape and normalization of the halo mass function. From 
Figure [12] the dependence of cosmology on this uncertainty 
is nearly negl i gible. In the analysis of maxBCG counts 
by iRozo et al.l (|2010), the dominant systematic uncertainties 
are miscentering and the normalization of the weak lensing 
mass calibration. Miscentering has less of an effect on M/N 
than on abundance analysis. The weak lensing measure- 
ments of ISheldon et al. (2009b) and halo mass estimates of 
iJohnston et al.l (120071) are used to calibrate the mass-richness 
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FIG. 16. — Left Panel: Comparison of the M/N results to parameter constraints from maxBCG cluster counts from|Rozo et al. (2010). Both results assume 
WMAP7 priors. It is important to note that cluster abundances are not used anywhere in our analysis. The systematic uncertainties that Rozo et al. are most 
sensitive to are the ones on which M/N have little dependence. Thus, the constraints from these two studies are roughly independent. Right Pane/:Comparison 
of this work to other cluster abundance studies. All contours are 95% and incorporate WMAP7 priors. The M/N technique is consistent and competitive with 
cluster abundance studies while using complementary information to cluster counts. 



relation in the lRozo et al.l (1201 Ol) results, thus the uncertainty 
in the weak-lensing calibration is correlated between the M/N 
and abundance analyses. However, this uncertainty is sub- 
dominant in the M/N analysis. Removal of this uncertainty 
entirely would only change the prior on / sys from 0.227 to 
0.225. A larg er error in the weak lensing calibration than that 
suggested by iRozo et all d2010l) would affect both measure- 
ments. A fully robust analysis of the maxBCG cluster catalog 
should simultaneously model NiNooo), w p (r p ), and M/N, but 
given the different dependences of the techniques on the com- 
mon systematic uncertainties, the two results can be consid- 
ered nearly independent. Combined constraints using M/N, 
cluster counts, and CMB data are listed in Table 5. For non- 
optical cluster counts, such as detections through X-ray or 
Sunyaev-Zeldovich observations, the dominant observational 
systematics are the absolute normalization and scatter of the 
mass-observable relation. Because a sat is close to unity for 
all galaxy samples, the M/N function is nearly independent 
of cluster richness and thus of scatter in mass versus richness. 
Figure QT| demonstrates the insensitivity of our results to our 
assumptions about the uncertainty in this scatter. 

If the theoretical uncertainties on e„, e&, and e/, were sig- 
nificantly smaller, our best fit value for as would be ~ 0.91. 
This would be in tension with the WMAP7 cosmological con- 
straints as well as cluster abundance results. This tension 
could be relieved if the true halo bi a s func tion had a higher 
normalization than the iTinker et al.l (1201 Ol) result by ~ 3%, 
or a scale-dependent bias that is closer to scale-independent 
relative to the non-linear dark matter clustering. 

6.2. Prospects for the Future 

Substantial improvements in the cosmological constraints 
presented in this paper can be achieved by reducing the sys- 
tematic uncertainties in the current analysis. The theoreti- 
cal uncertainties most applicable to our analysis, namely the 
large-scale bias and scale-dependent bias of dark matter ha- 
los, can be reduced with improved numerical simulations. Al- 



though state-of-the-art just two years ago, the N -body simu- 
lation set analyzed in ITinker et al.l (l2008al) and ITinker et al.1 
( 120101) has been replicated many times over (in volume and 
particle number) by ongoing numerical studies. The simula- 
tions of the Las Damas project 6 can currently achieve marked 
improvement in the halo mass function (McBride et al., in 
preparation) and the bias statistics for typical ACDM cos- 
mologies. A fully robust theoretical framework requires un- 
derstanding of halo populations in dark energy cosmologies, 
although we expect the impact of dark energy t o be small at 
fixed fl m , og, and power spectrum shape (e.g., iKuhlen et alj 
l2005h . 

For the systematic uncertainties in our measurements, there 
exists a simple solution to remove the uncertainty on the evo- 
lution of the luminosity function and the HOD between the 
clustering samples and the sample of clusters: measure the 
clustering from the same volume of space as the maxBCG 
sample. A spectroscopic sample of luminous red galaxies 
(LRGs) with n ? ~ 1 x 10" 4 extends to z = 0.36 (see, e.g., 
iKazinet alJfeOlOl) . Nearly every maxBCG cluster above a 
richness limit of A^oo = 9 (M sa 3 x 10 13 h~ l M Q ) should con- 
tain one or more LRGs (IZheng et al.ll2009l) . thus this analy- 
sis can be repeated on a sample of LRGs with an uncertainty 
on / sys than only comes from the weak lensing calibration. 
From the low number density of LRGs, as well as the halo 
occupation analysis of lZheng e t al. (2009), we expect (N^m 
for LRGs to be roughly 1 /4 that of the bright sample ana- 
lyzed here, increasing the statistical errors on (N sa t)M from 
~ 5% to 10%, but the size of the maxBCG cluster sample 
analyzed can b e increased by a factor of two relative to the 
sample used in ISheldon et alj d2009bl) and in this paper. In 
the near future, the Baryon O scillation Spectroscopic Survey 
(BOSS: ISchlegel et alJlBoH ) will increase the number den- 
sity of LRGs at z < 0.4 by a factor of three relative to the 
DR7 LRG number density. All of this points to luminous 

6 |http : //lss ■ phy . vanderbilt . edu/lasdamas/| 
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red galaxies as a logical and straightforward extension of the 
M/N analysis. 

In the longer term, there are numerous photometric surveys 
planned that will produce vast numbers of optically-detected 
galaxy clusters out to z > 1, including the Dark Energy Sur- 
vey, Pan-STARRS, and the Large Synoptic Survey Telescope. 
These surveys are also designed for weak lensing in order to 
obtain cluster masses as well as cosmic shear to measure the 
amplitude of dark matter clustering. These surveys will yield 
precise measurements of the angular clustering of L» galax- 
ies. The M/N technique can be applied to these data with lit- 
tle modification from the analysis presented here. Due to the 
complementary nature of the M/N approach to cluster abun- 
dances, this technique will enhance the measurements on the 
growth and expansion history of the universe and tighten con- 
straints on the equation of state of dark energy. For photomet- 
ric data, the dominant systematic will be uncertainties in the 
photometric redshift errors, and it will be important to quan- 
tify how the effect of photometric redshift uncertainty corre- 
lates M/N and cluster abundance data. As opposed to cluster 
counts, which are sensitive both to the growth of structure and 
to the expansion history of the universe, the M/N approach 
is primarily sensitive to growth because the volume element 
does not enter directly into the M/N predictions. Thus the 
combination of these two approaches may break the degener- 
acy between dark energy and models in which gravity is mod- 
ified at large scales in order to explain late-tim e acceleration 
(e.g jKnox et alJl2006tlHuterer & Linderj|2007l) . 

7. SUMMARY 

We have presented measurements of the mass-to-galaxy- 
number ratio (M/N) within maxBCG clusters in the SDSS. 
Combined with measurements of the projected clustering of 
galaxies, we have demonstrated that these data are a sensitive 
probe of the cosmological parameters il m and erg, while be- 
ing insensitive to other cosmological parameters that mainly 



enter into the shape of the linear matter power spectrum. Our 
theoretical analysis is based on the halo occupation distribu- 
tion (HOD). We have incorporated uncertainties in the the- 
oretical modeling, including uncertainties in the halo mass 
function, the halo bias function, and the scale-dependent bias 
of halos. We are most sensitive to uncertainties in the halo 
bias function and to the possibility of evolution between the 
mean redshift of the cluster sample, z ~ 0.25, and the mean 
redshift of the clustering sample, z ~ 0.1. After incorporat- 
ing all these uncertainties in our analysis, we find fij^ 5 as = 
0.465 ±0.026, with individual constraints of Cl m = 0.29 ±0.03 
and cr 8 =0.85 ±0.06. Combined with current CMB data, these 
constraints are tt m = 0.290 ±0.016 and cr 8 = 0.826 ±0.020. 
These constraints are consistent with and comparable to those 
obtained by cluster abundances, even though abundance in- 
formation is not used in this analysis. Future investigations 
could reduce the systematic uncertainties in this analysis and 
thereby tighten parameter constraints by a factor of two or 
more. 

The systematic uncertainties that are most important for 
cluster abundance studies, namely theoretical uncertainties in 
the halo mass function and scatter in the mass-observable re- 
lation, have negligible impact on the constraints achieved with 
M/N. Thus the combination of these two techniques can pro- 
vide a unique probe of the growth and expansion history of 
the universe from future photometric galaxy surveys. 
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APPENDIX 

A. MODIFIED SCALE-DEPENDENT BIAS 

In this paper we use a modified form of the scale-dependent bias function presented in lTinker et al.l (120051) . That function was 
calibrated on a series of N-bo dy simulations in which the halos were identified using the friends-of-friends percolation algorithm 
(FOF; e.g.. lDavis et al.|[T985l) . with a linking length of 0.2 times the me an interparticle separation . The halo mass function and 
halo bias function used here use the spherical overdensity algorithm (SO: lTinker et al.l2008all20Tob . In FOF, nearby halos can be 
linked together and labeled as a single object. In SO, this rarely happens. In the SO algorithm implemented by Tinker et. al., halos 
are allowed to overlap so long as the center of one halo is not within the virial radius of anoth er. Thus the s mall-s cale clustering 
of halos in these two scenarios should not be the same. For pair separations r > R\ + Ro, the iTinker et al.l (120051) function is an 
adequate description of the scale dependence of bias. For pairs in the overlap regime, r < R\ +R2, we find that the bias with 
respect to the non-linear matter distribution is nearly constant. Thus we adopt 



b 2 (M,r)-- 



b 2 (M) t' +L17 ^W] 

" ylV1> [l+0.696„(r)] 2l » 

h 2 (M \ [l+i.ne^Zfihaio)] 149 

° V V1 > [l+0.696„(2Rhui u )r- 09 



if r >= 2fl ha i 
if r < 2R h . dlo 



(Al) 



as our function al form. At r <_fthai , b(M, r) = by definition, but this is enforced in the halo exclusion of the two-halo term. See 
Appendi x B of rHnker et al.l (120051) for full details on the analytic model, with the minor differences described in ^3.11 

Figure IC17I shows the halo a utocorrelation functions for five mass bins. The data points are taken from two simulations 
described in Tinker eTail (120101) : L1000W and H384. Results are shown here for halos with N > 400 particles pe r halo. The 
curves are calculated assuming £(r) = b 2 (M, r)^ m (r), where £ m (r) is the non-linear matter correlation function given by lSmith et al.1 
J20031) . The dotted curves show the original function and the solid curves show the modification in eq. (lAlb . It is clear that eq. 
dAU is a better description of halo-halo clustering at small scales. The fit is not perfect, however, which is why we have adopted 
a 15% error on the deviation of the bias from scale independence. 
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B. EVOLUTION OF THE HOD IN SEMI-ANALYTIC MODELS 

One of our primary uncertainties is whether the HOD measured at z = 0.25 should be the same as that inferred from clustering 
at z = 0.1. In a separate paper we will demonstrate that the halo occupation of galaxies inferred through the subhalo abundance 
matching paradigm shows negligible evolution for samples defined at a fixed number density (Redd ick et al., in prepar ation). 
We have also investigated the evolution of the HOD in th e semi-analytic galax y formation model of iBower et al.l (|2006;i), built 
upon the high resolution Millennium N-body simulation dSpringel et al.ll2005l) . This model compares well with the observed 
luminosity an d color distribution of galaxies, and it does a reasonable job reproducing the dependence of galaxy properties on 
environment dBaldry et al.ll2006h . Figure |CT~81 shows the evolution of the HOD for galaxies in the SA model for the same three 
number density thresholds as our z = 0.1 clustering samples. In each panel, the upper frame shows the HOD at z = 0.25, z = 0.1 
and z = 0. The lower frame shows the ratio of (N)mz=o.25 to (N)m at z = 0.1 and z = 0. For the halo mass range probed by the 
maxBCG sample, and over the redshift baseline of our observations, (N)m varies by 5-10%, depending on the galaxy number 
density. We incorporate a 10% uncertainty in the evolution of the HOD in our analysis. 

C. TESTS WITH MOCK GALAXIES 

We have tested our methodology on mock galaxy distributions created on a high-resolution N-body simulation. The simulation 
is H384 from the simulation set analyzed in iTinker et al.l (l2008al 120101) . The simulation consists of 1024 3 particles evolved in 
a volume 384 h~ l Mpc per side. The mass resolution of this simulation is high enough to resolve halos well below M m ;„ of the 
01 M r < -19.5 galaxy sample. The simulation cosmology is (il m ,a^,ilb,h(),n s ) = (0.3,0.90.04,0.7, 1.0). The halo catalogs are the 
same as those in the Tinker et. al. mass function and bias relation analyses, using the SO halo finder. 

Our procedure for creating the mock data was as follows. First, we fit the w p (r p ) data alone assuming the same cosmology 
as the simulation. The halos within the simulation are populated with galaxies according to the best-fit HODs. For each halo, 
a maxBCG richness A^oo is assigned according to the mass-richness relation given in eq. ( TTOb (including scatter). The halos are 
stacked in the same richness bins given in Table 2, and the mean value of M/N for each bin is calculated directly from the set of 
halos in each bin; i.e., we do not create shear profiles or projected number density profiles from the galaxy distribution. Our test 
was directed at the methodology of combining w p (r p ) and M/N itself, rather than the techniques for measuring these quantities, 
which are tested elsewhere and are included in the systematic error budget in our analysis. We include all the same systematic 
uncertainties in this analysis — halo statistics, mass-richness relation, cosmological parameters — except for / sys , which we fixed 
at unity and do not allow to vary. We obtain the error bars on w p (r p ) from performing the jackknife technique on the simulation 
itself, but for M/N we adopt the same (proportional) errors as on the maxBCG data. 

Figure |CT9l shows the cosmological constraints produced by the mock analysis. The error contours are centered on the in- 
put cosmology. The banana curve is similar to the results from the actual data but with a stronger constraint on the cluster 
normalization due to the tight prior on / sys . 
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FIG. C17. — Halo-halo correlation functions from two high reso luti on N-body simulat ions (the L1000W and H384 simulations utilized to calibrate both 
the halo mass function and halo bias relation in Tinker et al. 2008a and Tinker et al. 2010). The three low-mass bins are from the H384 simulation and the 
three high-mass bins are from the L 1 000W simulation. The abrupt increase in c lustering amplitude between the two simulation results is due to the change in 
cosmological model. The dotted line shows the scale -dependent bias fun ction of Tinker et al. (2005). The solid curves show our modification to this formula to 
account for the change in halo definition between the Tinker et al. 1 2005) simulations and the more recent calibrations. See text for details. 
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FIG. C18.— Evolution of the HOD in the lBower etaH f200f3) semi-analytic galaxy formation model. The right, middle, and left panels show results for galaxy 
samples that match the number densities of our M r < -21, o i M r < —20.5, and o i M r < —19.5 samples, respectively. The lower frame in each panel shows the 
ratio with respect to the z = 0. 1 HOD. 
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FIG. C19. — Cosmological constraints on m ock data created from an N-body simulation 384 h Mpc per side (the H384 simulation utilized to calibrate both 
the halo mass function and halo bias relation in Tinker et al. 2008a and Tinker et al. 2010). The mock analysis incorporates uncertainties on e, u t r , and ej, but 
assume that / sys = 1 with no error. The circle indicates the input cosmology. 
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